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Introduction: Circulation Control (CC) airfoil 

~ For the next generation of passenger aircraft ~ 


Circulation control by “Coanda jet” 

Jet flow attaches on a rounded trailing 
edge of an airfoil 

Circulation is increased 
Lift is enhanced 



Hydrogen bubble flow visualization of 
a CC airfoil (NASA LaRC, 2002) 


CC for Aircraft applications 

- Lift enhancement (by single jet-blowing) 

- Maneuver support (by dual jet-blowing) 


Coanda effect: 

The tendency of a fluid jet to stay attached to an Concept image of hybrid-wing-body 

adjacent curved wall, named after Henri Coanda aircraft employing CC devices 
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Introduction: Separation of turbulent Coanda flow 

- A major difficulty in current RANS Simulations ~ 


Coanda jet over a cylinder Coanda jet over a rounded trailing 
(Exp. by Wygnanski et al.) Edge of a CC airfoil 



No external flow 


Transitional 

. Gortler) 


Jet flow characteristics are sensitive 
to the transition process 



Turbulent boundary layer 


Immediate transition 
to turbulence 


Turbulent boundary layer * 

Jet flow develops to fully turbulent at the jet exit 


Jet plenum 


Coanda sitrfac 
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Objectives of the study 


1 . To investigate detailed physics (flow structures and statistics) of the 
fully turbulent Coanda jet applied to a CC airfoil, by using LES 

2. To compare LES and RANS results to figure out how to improve the 
performance of existing RANS models for this type of flow 
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Airfoil configuration 



c = 218.3mm 


- 20%-thick, non-cambered, elliptic-leading-edge CC airfoil 

(designed by Dr. Englar’s research group at GTRI - currently tested at NASA Langley) 

- Two independent jet plenums for the upper and lower sides 

(lower jet slot is closed in this study - may be used for “dual blowing” in future studies) 

- Chord Reynolds number: 0.49 million 

- Wavy transition strip attached near the leading edge 
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Numerical methods (LES) 


Incompressible Navier-Stokes solver “CDP” 

(developed at the Center for Turbulence Research, Stanford University) 

- Unstructured, finite-volume solver 

- Energy-conservative, 2nd-order central difference scheme 

- Fully-implicit, 2nd-order time integration scheme 

- Dynamic Smagorinsky model for the subgrid-scale (SGS) stresses 

- Running on massively parallel supercomputers (256 CPU’s used) 
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Computational domain & Boundary conditions 


Airfoil chord 
c = 218.3 mm 

(Re c = 0.49 x 10 6 ) 

4.651c 


Wind tunnel geometry: NASA Langley Basic Aerodynamics Research Tunnel (BART) 

Spanwise (periodic) domain size: 14 mm = 0.0641c = 27.8 h 

- The domain wide enough to study turbulent structures in the Coanda jet 

- 3D RANS study has shown little sidewall effects at low jet-blowing case 



i 

k Top wall (inviscid) 

Inlet 


2.326c Outlet 

(14= 34 m/s) ^ 


(convective) 

: — ► 


c 

2.326c 



r Bottom wall (inviscid) 

< ► 

< ► 


3.619c 10c 
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Jet blowing conditions (at the jet exit) 


Jet- blowing rate 


Low 

High 


Uj [m/s 


105.3 

173.4 



max 



« 135 

rs 216 


Re j 

prf 4500 
« 7200 


UjA 

0.044 

0.120 


Uj : Bulk jet velocity (mean of the time-averaged velocity profile) 

Uj max : Maximum jet velocity (maximum of the time-averaged velocity profile) 
Rej : Jet Reynolds number ( = U jmax h/ v) 

C M : Jet momentum coefficient ( = 2 U?h / Uj c ) 



c - 21 8.3mm 
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Mean flow profiles (at the jet exit) 


Streamwise velocity Reynolds stresses 



Jet oscillating due to alternating 
(von Karman type) vortex shedding 
behind the thin jet blade 

0,91 0^2 0^3 

x/c 
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Computational grids 


y/c 0 


-0.2 

I 

# 

-0.4 

/ 

— 

/ 



0 0.2 
Wavy transition strip 


Multi-block 
structured grids 
for the whole domain 

3 levels of grids for 
the Coanda region 
(N0 x Nr x Nz) 

Coarse 

400 x 160 x 128 
Medium 

800 x 160 x 128 
Fine 

800 x 160x256 


Total number of grid points: 116 million (Fine grid) 
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Wavy transition strip 


Z 


< 




Wavy strip directly meshed 
using multi-block structured grids 

Strip height: 0.56 mm 

Strip width: 3.0 mm 

Wavelength: 7.0 mm 

Amplitude: 2.0 mm 
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Grid resolution (Fine grid) 


Coanda surface Resolution [nun] Resolution in wall units 

(0 = 10° ~ 170°) (0 = 15°) (0 = 30°) (0 = 45°) 


Chord wise 

0.0902 

Pi 95 ft! 65 

Pi 30 

Wall- normal 

0.00127 

ps 1.4 ft; 1.0 

Pi 0.5 

Spanwise 

0.0547 

pi 55 Pi 40 

Pi 18 

Airfoil surface 

Resolution [mm] 

( xjc — 0.5) 

Resolution in wall units 

(upper, xjc = 0.5) (lower, xjc = 0.5) 

Chordwise 

0.4906 

pi 70 

S3 45 

Wall- normal 

0.00254 

Pi 0.4 

« 0.25 

Spanwise 

0.0547 

Pi 7.5 

ss 5.0 


SGS eddy viscosity: Up to about 5 times larger than the molecular 

viscosity ju (just downstream of the jet exit, 
where RANS eddy viscosity is about 40 to 50 ju) 
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LES results: Flow around the airfoil 


Instantaneous 02 
streamwise 0 1 
velocity 

y/c o 
-o.i 
- 0.2 



Time-averaged 

pressure 

coefficient 
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LES results: Flow around the airfoil 


Instantaneous spanwise vorticity 



x/c 
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LES results: Flow around the airfoil 





Skin friction coefficient c, 

J 0 002 0.04 0.06 008 
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LES results: Self-sustained transition 


Leading edge 
if 



I L 


0 0.1 

x/c 


Turbulence was sustained with 
no “inlet disturbances” given in 
the present LES (disturbances 
were given to the whole domain 
only at the initial stage) 

Upper side 

Spanwise modulation 

I 

Transition to turbulence 
sustained around x/c = 0.1 

Lower side 

3D separation behind the strip 

I 

Transition to turbulence 
sustained around x/c = 0.4 
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LES results: Flow around the jet exit 
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LES results: Coanda jet profiles 


Top: Streamwise velocity 
Bottom: Reynolds shear stress 
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Angle 

(degrees) 



LES results: Vortical flow 
structures in the Coanda jet 

Isosurfaces of Q 

(2nd Invariant of velocity gradient tensor) 
colored based on velocity magnitude 
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LES results: Backward-tilted hairpin vortices 



between the legs outside the hairpin 

- Backward-tilted (i.e., head of each hairpin is located upstream of its legs) 

- Located above the high-momentum jet flow 

- Creating a strong upwash between the legs 

-> Lifting the high-momentum flow upward -> Turbulent mixing enhanced 
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LES results: Plots of velocity fluctuations 


0 = 30 degrees 
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LES results: Plots of velocity fluctuations 


0 = 30 degrees 




Uq/Uq 


m NASA 

WASA Postdoctoral Program 

administered by Oak Ridge Associated Universities 


Page 23/27 


Comparisons between LES and RANS 


Mean pressure distributions 
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Comparisons between LES and RANS 


RANS (S-A) 

Cp 

-4 -3 *2 *1 0 


RANS (SST) 

Cp 

-4 -3 - 2-10 


LES 


Cp 

-4 -3 *2 *1 0 



0.1 

0.05 

y/c o 
- 0.05 
- 0.1 




Jet separation: 75.0 deg. 
Bubble size: 0.058c 

Lift coefficient: 1 .85 


69.0 deg. 
0.080c 
1.60 


69.5 deg. 
0.060c 
1.36 
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Comparisons between LES and RANS 


Streamwise 

velocity 


Reynolds 
shear stress 
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Conclusions 


High-resolution LES of a turbulent Coanda jet (applied to a circulation 

control airfoil) was performed and was compared with RANS results 

LES results: 

1 . Pressure distributions agreed well with the preliminary experiments 

2. Many “backward-tilted” hairpin vortices were observed in the outer 
shear layer of the jet; the hairpins lift high-momentum flow upward 

Comparisons between LES and RANS: 

3. S-A and SST models predicted a larger circulation and a higher lift, 
even though SST model predicted a correct jet separation location 

4. Both models predicted a smaller jet spreading rate than the LES as 
the eddy viscosity was too small in the outer shear layer of the jet 
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(backup slide) 



Similar but opposite to “forward-tilted” hairpins in boundary-layer flows 

Boundary-layer flows: dil/dy > 0 

Wall-jet flows: dU/dy < 0 
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